Pinpointing Cosmic Ray Propagation With The AMS-02 Experiment 
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The Alpha Magnetic Spectrometer (AMS-02), which is scheduled to be deployed onboard the In- 
ternational Space Station later this year, will be capable of measuring the composition and spectra 
of GeV-TeV cosmic rays with unprecedented precision. In this paper, we study how the projected 
measurements from AMS-02 of stable secondary-to-primary and unstable ratios (such as boron-to- 
carbon and beryllium-10-to-beryllium-9) can constrain the models used to describe the propagation 
of cosmic rays throughout the Milky Way. We find that within the context of fairly simple propaga- 
tion models, all of the model parameters can be determined with high precision from the projected 
AMS-02 data. Such measurements are less constraining in more complex scenarios, however, which 
allow for departures from a power-law form for the diffusion coefficient, for example, or for inhomo- 
geneity or stochasticity in the distribution and chemical abundances of cosmic ray sources. 



I. INTRODUCTION 

Despite nearly a century of observational and theoret- 
ical progress, the origin of the cosmic ray spectrum re- 
mains a major puzzle of modern astrophysics. The task 
of identifying the sources of these particles is complicated 
by the non-trivial processes involved in cosmic ray prop- 
agation. Although cosmic ray composition and spectrum 
measurements have taught us a great deal about the ac- 
celeration and propagation of cosmic rays, we still lack a 
detailed and self-consistent understanding of how these 
particles are produced, and how they travel through and 
interact with the interstellar medium. 

Below the spectral feature known as the knee (E ~ 
10 15 eV), the bulk of the cosmic ray spectrum is believed 
to be of galactic origin. Non-rclativistic shocks occurring 
in supernova remnants seem to be likely sources, and 
are predicted to accelerate cosmic rays with a power-law 
injection spectrum Q oc E^ 1 with 7 ~ 2 [1]. At higher 
energies, even less is known about the origin of the cosmic 
ray spectrum. In this work, we focus solely on galactic 
cosmic rays at energies well below the knee. 

Our understanding of how cosmic rays propagate 
through the Milky Way is informed largely by mea- 
surements of the spectra of various cosmic ray species 
as observed at Earth. In particular, by comparing the 
spectrum of particles produced in cosmic ray accelera- 
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tors (primaries) to those that are produced by inelastic 
processes during propagation of primary particles (sec- 
ondaries), we can learn about the mechanisms involved 
in cosmic ray propagation. While stable secondary-to- 
primary ratios (such as boron-to-carbon and antiproton- 
to-proton) provide information that can be used to con- 
strain the effective column density cosmic rays pass 
through before reaching the Solar System, unstable ra- 
tios (such as beryllium-10-to-beryllium-9) are useful in 
constraining the time interval since spallation. Combina- 
tions of such observations make it possible to constrain 
the basic properties of relatively simple cosmic ray prop- 
agation models. 



To date, some of the most precise cosmic ray mea- 
surements over the GeV-TeV energy range have been 
made by the CREAM (boron-to-carbon) [2], PAMELA 
(antiproton-to-proton) [3], ISOMAX ( 10 Be-to- 9 Be) [4], 
and HEAO-3 (Subiron-to-iron, boron-to-carbon) [5] ex- 
periments. These measurements have been used to place 
fairly stringent constraints on the parameters of the un- 
derlying cosmic ray propagation model [6, 7]. In this 
article, we extend this approach to include data antici- 
pated from the Alpha Magnetic Spectrometer (AMS-02) 
experiment, which is scheduled to be deployed on the 
International Space Station in 2010. With its greater ac- 
ceptance and superior particle identification relative to 
previous experiments, measurements from AMS-02 are 
expected to dramatically improve our understanding of 
the processes involved in galactic cosmic ray propagation. 
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II. COSMIC RAY PROPAGATION IN THE 
MILKY WAY 

Once injected from their sources into the interstellar 
medium, charged cosmic rays — unlike photons or neu- 
trinos — undergo a number of processes potentially ca- 
pable of significantly altering their spectra (for a recent 
review, see Ref. [8]). The Galactic Magnetic Field, in 
particular, is responsible for deflecting charged particles, 
leading them to diffuse gradually throughout the Galaxy, 
following paths resembling a random walk. Particles with 
greater energy, and therefore rigidity, diffuse more ef- 
ficiently and tend to escape the Galaxy more quickly, 
whereas less energetic cosmic rays are typically confined 
by the Galactic Magnetic Field for a greater duration. 

An essentially inevitable consequence of high energy 
particle scattering in the turbulent magnetic field is 
stochastic acceleration, also known as diffusive reaccel- 
cration [9] . This mechanism gives rise to diffusion in mo- 
mentum space with a diffusion coefficient determined by 
the spatial diffusion coefficient and the Alfvcn velocity, 
which represents the typical velocity at which magnetic 
irregularities propagate in the interstellar medium. 

Other potentially important effects to consider include 
galactic winds, which may result in the convection of par- 
ticles away from the Galactic Plane, as well as various 
energy loss processes. Such energy losses occur as a re- 
sult of the cosmic rays traversing the galactic medium, 
which is permeated with gas, radiation fields, and mag- 
netic fields. In the case of nuclei Coulomb and ion- 
ization energy losses exist, but they play only a minor 
role in their propagation. On the other hand, GeV elec- 
trons and positrons lose significant quantities of energy 
through inverse Compton and synchrotron processes; at 
lower energies, ionization, Coulomb interactions, and 
bremsstrahlung processes may also be relevant. Further- 
more, the decays of unstable, radioactive species must 
be taken into account, including the introduction of any 
relevant decay products. Lastly, cosmic ray spallation on 
the interstellar medium can lead to the extinction of the 
incident particle and to the creation of a secondary flux 
consisting of gamma rays from neutral pion decay, elec- 
trons, positrons, protons, antiprotons and nuclei. The 
rate at which spallation occurs is fixed by the nuclear 
cross sections involved, and by the distribution of gas 
(mainly H and He) present in the Milky Way. 

The transport equation that describes all of the above- 
mentioned processes for a cosmic ray species i (with 
atomic number Zi and mass number Ai) is given by [1, 8] 
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sity of particles of species i per unit momentum, R = 
pel (-Zi|e|) is the rigidity, D xx (x, R) is the spatial diffusion 
coefficient, V c (x) is the convection velocity, D pp (x,R) is 
the momentum space diffusion coefficient, Td.i is the de- 
cay time, T sPi i is the spallation time, and — p\{x,p) is the 
energy loss rate. Notice that the source term Qtot,i in- 
cludes the primary injection spectrum Q in j cx i?~ 7 (typ- 
ically assumed to be from acceleration in supernova rem- 
nants) as well as the products of decay and spallation 
of heavier cosmic ray species. Qtot,i may also include 
exotic primary contributions produced, for example, by 
dark matter annihilations taking place in the halo of the 
Milky Way. Such components are not considered in this 
study, however. 

The standard approach used to solve Eq. (1) is to as- 
sume that the cosmic rays are in steady state, such that 
drii/dt = (or to solve iteratively the full equation un- 
til an approximate steady state is reached). To further 
simplify the problem, it is common to adopt a cylindrical 
diffusive region (of radius r max and half-thickness L) in- 
side which the diffusion coefficient does not vary with lo- 
cation, and such that the number density of cosmic rays 
approaches zero at the boundaries. The half-thickness 
of the cylindrical diffusive region, L, is generally much 
larger than the half-thickness of the galactic disk h ~ 0.1 
kpc. 

To proceed, one may either apply (semi-)analytical 
methods (e.g. [10-12]) or make use of numerical codes 
(e.g. GALPROP [13, 14] or DRAGON [7, 15]). In the 
present work, we use GALPROP v50.1p [13] which solves 
Eq. (1) by assuming a homogeneous power-law coeffi- 
cient, 

D XX {R) = (v/c)D 0xx (R/R ) a , 

in a cylindrical diffusive region as described above where 
the Sun sits at x = (r ,z ) = (8.5,0) kpc. GALPROP 
implements a realistic interstellar hydrogen distribution 
and interstellar radiation field based on state-of-the-art 
surveys, and makes use of an up-to-date nuclear reac- 
tion network. The primary sources of cosmic rays are 
assumed to have a unique single or double power-law in- 
jection spectrum and isotope composition. Moreover, the 
distribution of sources can be specified but is typically 
chosen such that the EGRET gamma ray data are well 
reproduced. Convection, if implemented, is assumed to 
move particles along the vertical direction, perpendicular 
to and away from the disk and with a linear profile, 

V c (x) = sgn(z)(V c . + \z\dV c /dz)e z . 

Diffusive reacceleration with a given Alfven velocity, va, 
is also included, and is modelled with a diffusion coeffi- 
cient 



D pp (R) 
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The local interstellar flux of species i, $j, follows from 
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the solution of Eq. (1) through 

cA 

$i(x Q ,T) = -^m(x Q ,p), 

where T is the kinetic energy per nucleon. Finally, this 
flux is modulated by the solar wind before arriving at 
the top of the atmosphere, as calculated using the force 
field approximation [16]. Further details pertaining to 
the GALPROP package are described in Refs. [13, 14]. 

III. PROSPECTS FOR AMS-02 

The second and final version of the Alpha Magnetic 
Spectrometer (AMS-02) [17, 18] is a large acceptance 
cosmic ray detector scheduled to be placed onboard the 
International Space Station in 2010. Over its mission 
duration of at least three years, it will measure with un- 
precedented statistics and precision the spectrum of cos- 
mic rays over an energy range of approximately 100 MeV 
to 1 TeV [19, 20]. The AMS-02 instrument will be able to 
detect and identify nuclei as heavy as iron (Z < 26) with 
rigidity up to 4 TV [18], and separate isotopes of light 
elements (namely, H, He and Be) over a kinetic energy 
range of 0.5—10 GeV/n [21-23]. High precision measure- 
ments of the ratios B/C and sub-Fe/Fe (D/p, 3 He/ 4 He 
and 10 Be/ 9 Be) up to energies of-1 TeV/n (-10 GeV/n) 
are anticipated. Due to a high level of proton rejection, 
positron and antiproton spectra will also be measured 
with unprecedented precision [22, 24]. 

The principal goals of the AMS-02 experiment [18, 25] 
include searches for primordial anti-matter among cos- 
mic ray nuclei and signatures of dark matter annihila- 
tions or decays. Additionally, however, AMS-02 will con- 
tribute considerably to our understanding of the origin, 
acceleration, and propagation of galactic cosmic rays. In 
fact, being the largest acceptance (— 0.45 m 2 sr) space- 
based magnetic spectrometer, AMS-02 will bypass the at- 
mospheric systematics which affect baloon-borne exper- 
iments and simultaneously feature an acceptance more 
than two orders of magnitude above the — 0.002 m 2 sr 
of the PAMELA satellite [26]. Besides largely improved 
statistics, AMS-02 will also provide cosmic ray flux mea- 
surements up to — TeV/n and separate Be isotopes up to 
— 10 GeV/n while PAMELA can only reach a few hun- 
dred GeV/n and separate H and He isotopes. Therefore, 
a precise B/C ratio over a wide energy range and high 
energy 10 Be/ 9 Be measurements will be obtained. These 
data are extremely useful in constraining cosmic ray dif- 
fusion parameters. Of particular importante will be the 
ratio 10 Be/ 9 Be since as of today there are no data points 
above 3 GeV/n. 

Focusing on the ratios B/C, 10 Be/ 9 Be and p/p, we can 
estimate the systematic and statistical errors of AMS- 
02. To begin, the rigidity resolution of AMS-02 will 
be AR/R -1 - 2% at -10 GV, for both protons and 
He nuclei, and around 20% at -0.5 TV (-1.0 TV) for 
protons (He) [22, 27, 28]. Assuming similar capabilities 



for heavier nuclei as well, we take a conservative value 
AR/R — 20%, which in the case of relativistic parti- 
cles translates directly into a kinetic energy resolution 
AT /T — 20%. Such resolution allows logarithmic bins of 
width log 10 — 0-087, or 11 — 12 bins per decade. 

In the following we assume 10 bins per decade of kinetic 
energy regardless of the cosmic ray species. 

In order to compute the statistical errors associated 
with the ratio Ni/Nj, we need the number of i and 
j particles detected, Ni = Ci acci <E>j ATAti (and like- 
wise for j), Ci being the efficiency, acci the geometri- 
cal acceptance of the instrument, and Ati the operating 
time. Then, MN ffl£ 3tat = l/y/Ni + 1///V-. The ge- 
ometrical acceptance is a function of the particle type; 
we adopt acc B = accc = acc Be = 0.45 m 2 sr [18, 29], 
acc p = 0.3 m 2 sr [23], and accp — 0.160(0.033) m 2 sr for p 
momenta 1 - 16 (16 - 300) GeV [19, 29, 30]. Following 
Ref. [28], we fix e B = ec = 95% and all other efficiencies 
to 100%. Lastly, we consider one year of operation. 

As for systematics, we also estimate the errors associ- 
ated with the mismeasurement of the atomic number of 
cosmic ray nuclei. Using the full capabilities of the AMS- 
02 silicon tracker, the author of Ref. [28] used Monte 
Carlo simulations to estimate the level of misidentifica- 
tions, finding fewer than one percent for 2 < Z < 11. 

Conservatively, we take (AA ^"" < = (AA ^""' = 1%. 
The 10 Be/ 9 Be measurement, on the other hand, is more 
delicate since it relics not only on charge, but also mass 
separation. From Ref. [31], we take a mean mass res- 
olution for Be of Am/m — 2.5%. Requiring a sepa- 
ration of consecutive isotopes within 0.5 atomic mass 
units, this mass resolution results in misidentification of 
9 Be as 10 Be and vice-versa less often than / — 2.275%. 
Hence we use a systematic error for Nw Be given by 
{ANw Be ) syst /Nw Be = | - / + / N9 Be /Nio Be \, and anal- 
ogously for Ng Be . Clearly, the systematics become un- 
acceptable when either / Ng Be /Nio Be or / Nia Be /N9 Be 
approaches unity. Finally, for the antiproton-to-proton 
ratio, the dominant fraction of the systematic error comes 
from the p measurement. In order to confidently iden- 
tify antiprotons, the large background of protons and 
electrons must be rejected with high efficiency. In the 
multi-GeV energy range, at which p/p — 10 4 (e~ /p — 
10 2 ) [3, 30], the rejection power expected for AMS-02 is 
p : p - 10 5 - 10 6 (e- : p - 10 3 - 10 4 ) [29, 30], leading to 
a systematic error of order 10 19 1Q6 = 1 — 10% ( 10 3 1 ° 10 4 = 
1-10%). Consequently, we take (AA ^ 3 " 3f = 5%. We add 
all of the systematic errors discussed in this paragraph 
in quadrature to the statistical uncertainties. 

At this point we can forecast the AMS-02 measure- 
ments of various cosmic ray spectra. As a benchmark 
model, we adopt the best-fit propagation parameters 
found in Ref. [6]: D 0xx = 6.04- 10 28 cm 2 /s (at a reference 
rigidity of i? = 4 GV), L = 5 kpc, a = 0.41, v A = 36 
km/s, and no significant convection. The remaining spec- 
ifications are as in galdef_50p_599278 file [13], including 
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B/C for model (D_0xx=6.04x10 zs cm 2 /s, L=5 kpc, ct=0.41) 
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FIG. 1: Projections for ability of the AMS-02 experiment to 
measure selected stable secondary-to-primary and unstable 
cosmic ray ratios: B/C, 10 Be/ 9 Be and p/p. As a benchmark, 
we have considered the best-fit model of Ref . [6] . Systematic 
and statistical uncertainties are included in the AMS-02 error 
bars, and 1 year of data taking is assumed. The thick verti- 
cal line indicates the energy cut, T > 5 GeV/n, imposed to 
reduce the impact of solar modulation on our results. In the 
lower frame, the PAMELA measurement of the antiproton- 
to-proton ratio [3] is also shown. 



a distribution of cosmic ray sources optimized to meet 
EGRET gamma ray data, and a double power-law injec- 
tion spectrum — with indices 71 = 1.82 and 72 = 2.36 
below and above Rq — 9 GV — to reproduce low-energy 
cosmic ray data. The ratios B/C, 10 Be/ 9 Be and p/p cor- 
responding to this model (here after the true model) and 
calculated with GALPROP v50.1p are plotted in Fig. 1 
along with the projected error bars of the AMS-02 in- 
strument (including both systematic and statistical un- 
certainties). For detailed simulations of the capabilities 
of AMS-02, we refer the reader to Ref. [21, 31] and refer- 
ences therein. 

As shown in the lower frame of Fig. 1, the antiproton- 
to-proton ratio for the true model is somewhat lower than 
the values measured by the PAMELA collaboration [3], 
especially at low energies. Solar modulation, however, 
may have a significant impact on this ratio at such ener- 
gies. As the impact of solar modulation varies with re- 
spect to the time period observed, the PAMELA antipro- 
ton measurement is not necessarily expected to mimic 
that to be measured by AMS-02. To reduce the depen- 
dence on this effect, we apply an energy cut T > 5 GeV/n 
throughout our analysis. In our calculations, we have 
modulated the cosmic ray spectra with 4>p = 450 MV. 
Common values of the modulation parameter 4>p range 
from a few hundred MV up to over a GV (e.g. [32]); how- 
ever, 450 MV is a reasonable value for data taken around 
a solar minimum which presumably will be the case of 
AMS-02 first year. In any case, we stress that our results 
do not depend much on the solar modulation parameter 
since only energies above 5 GeV/n are considered. In 
the remainder of the work, unless otherwise stated, we 
shall use the modulated data set presented in Fig. 1 to 
perform our analysis. 



IV. CONSTRAINING PROPAGATION MODELS 

In this section, we attempt to estimate how well the 
projected AMS-02 measurements described in the pre- 
vious section will be able to constrain the propagation 
model parameter space. For the moment, we fix the 
Alfven speed to va = 36 km/s, neglect the effects of 
convection (V c ,o = dV c /dz — 0), and proceed in a fashion 
similar to Ref. [6] to run GALPROP 245 times, in a 7x7x5 
grid of the parameters (Dq xx , L, a) over the following 
ranges: D Qxx = 4.54-8.03- 10 28 cm 2 /s, L = 3.5-6.5 kpc 
and a = 0.39 — 0.43. Linearly-(Logarithmically-)spaced 
gridpoints were implemented for L, a (Dq xx ). An infill of 
3 points between consecutive gridpoints (corresponding 
to a reduction of the spacing by a factor 4) was performed 
and the relevant cosmic ray ratios for each additional 
propagation model were obtained through 3-dimensional 
interpolation of the GALPROP runs. The extended grid 
includes 25x25x17=10,625 different parameter sets. For 
each set, we calculate the % 2 using the projected B/C 
and 10 Be/ 9 Be AMS-02 measurements. 

Using the above described parameter scan and the pro- 
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FIG. 2: Constraints on cosmic ray propagation parameters 
from the projected AMS-02 B/C and 10 Be/ 9 Be measurements 
presented in Fig. 1. The solid lines delimit 1, 2, and 3 a 
regions using the B/C data set, whereas dashed lines refer to 
10 Be/ 9 Be. Here, the propagation parameters were varied in 
the ranges D 0xx = 4.54 - 8.03 ■ 10 28 cm 2 /s, L = 3.5 - 6.5 
kpc and a = 0.39 — 0.43. We have assumed va = 36 km/s, 
V c ,o — dV c /dz — 0, and have marginalized over a. 



jected B/C and 10 Be/ 9 Be presented in Fig. 1, we show in 
Fig. 2 the resulting 1, 2 and 3a regions in the L vs. D 0xx 
plane where we have marginalized over a. From this 
figure, one can immediately identify the complementar- 
ity between stable secondary-to-primary and unstable ra- 
tio measurements. Whereas stable secondary-to-primary 
ratios provide an approximate measure of the quantity 
L/Dq xx , unstable ratios help to determine L 2 /Do xx (for a 
fixed value of a). The combination of B/C and 10 Be/ 9 Be 
measurements can thus provide a determination of both 
L and Dq xx . Although the projected p/p data, shown 
in the lower frame of Fig. 1, introduces some additional 
information into the analysis, it provides a constraint re- 
gion with a similar shape but broader than that provided 
by the B/C data. For this reason, we do not include p/p 
in our chi-squares. 

Fig. 3 shows, for the same scan of propagation param- 
eters as used in Fig. 2, the 1, 2 and 3a contours from 
the combination of B/C and 10 Be/ 9 Be projected mea- 
surements presented in Fig. 1. In each frame we have 
marginalized over the parameter not shown. As this fig- 
ure demonstrates, the projected AMS-02 measurements 
of B/C and 10 Be/ 9 Be are sufficient (within the context 
of the simple models presently being considered) to de- 
termine the underlying propagation parameters with an 
accuracy of AD 0xx ~ 1.4 • 10 28 cm 2 /s, AL — 1.0 kpc, 
and Aa ~ 0.02 (at la). This precision is much greater 
than that obtained with present (pre- AMS-02) data; see 
Refs. [6, 7]. In particular the degeneracy between Dq xx 
and L is broken as can be seen in the upper frame of 
Fig. 3 where we have overplotted in dashed the 3a con- 
tour from Ref. [6]. 



o 

Q. 



10 



i — i — i — | — i — i — i — | — i — i — i — | — i — ,— 

+ true model 
X best-fit model 

- - - Simet & Hooper [6], pre-AMS02 (3a) 




28.2 



28.4 



28.6 



28.8 



29 



log10(D_0xx cms) 



Ax'fB/C+WBe) 



So.43 
0.425 
0.42 
0.415 
0.41 
0.405 
0.4 
0.395 
0.39 



-i— i— i— ,— i— i— r— i— r 
+ true model 
X best-fit model 




28.7 28.75 28.8 28.85 28.9 

log10(D_0xx cm' 2 s) 




6 6.5 

L [kpc] 

FIG. 3: Regions consistent (within 1, 2 and 3cr) with pro- 
jected B/C and 10 Be/ 9 Be AMS-02 data from Fig. 1 in the 
L vs. Dqxx, ot vs. Dqxx, and a vs. L planes. Here, the 
propagation parameters were varied in the ranges Do xx = 
4.54 - 8.03- 10 28 cm 2 /s, L = 3.5-6.5 kpc and a = 0.39-0.43. 
We have assumed va = 36 km/s, V c ,o = dV c /dz — 0, and have 
marginalized in each frame over the parameter not shown. In 
the top frame we show in dashed the 3<r contour from Ref. [6] . 
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(6.04, 5.000, 0.4100) 
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(6.04, 8.000, 0.4850) 
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(6.04, 8.997, 0.5000) 


2.17 





10 


(3.89, 5.623, 0.5000) 


12.78 



TABLE I: Best-fit models and reduced chi-squares for the 
combined B/C+ 10 Be/ 9 Be projected AMS-02 data set, for var- 
ious combinations of reacceleration and convection parame- 
ters. In each case, the data set used is that presented in 
Fig. 1. Propagation parameters were varied in the ranges 
D 0xx = 0.57 - 19.5 ■ 10 28 cm 2 /s, L = 1.22 - 20.48 kpc and 
a = 0.32 — 0.50. va and dV c /dz are given in units of km/s 
and km/s/kpc, respectively, and V c ,o = 0. The presence of 
significant convection will be highly observable to AMS-02. 
To a lesser extent, the value of the Alfven velocity will also 
be testable. 



Thus far, we have not considered any effects of con- 
vection and/or variations in the Alfven velocity from 
our default value of va = 36 km/s. As the quan- 
tity and quality of cosmic ray data improves, however, 
it will become increasingly possible to test these as- 
sumptions, and determine the related parameters. In 
order to study this possibility, we have repeated the 
procedure described in the previous paragraphs using 
the following values for the Alfven velocity and the 
convection velocity: (^[km/s], dV c /dz\km/s/kpc\) — 
{(36, 0), (15, 0), (0, 0), (0, 10)}, while leaving V c , = 0. We 
have focused here on Alfven velocities smaller than the 
reference value 36 km/s since some studies (e.g. Ref. [7]) 
find that lower values of va are preferred. For each com- 
bination of va and dV c /dz, we ran 343 GALPROP jobs in 
a 7x7x7 grid with ranges D 0xx = 0.57- 19.5- 10 28 cm 2 /s, 
L = 1.22 - 20.48 kpc and a = 0.32 - 0.50. Linearly- 
(Logarithmically-)spaced gridpoints were implemented 
for a (D 0xx , L) . An infill of 3 points between consecutive 
gridpoints was applied resulting in 25x25x25=15,625 dif- 
ferent models. Using the projected B/C and 10 Be/ 9 Be 
data sets shown in Fig. 1, we found the best-fit models 
given in Table I. As greater departures from our default 
assumptions are considered, the fits to the projected data 
become considerably worse. In particular, even modest 
(~10 km/s/kpc) amounts of convection lead to very poor 
fits to the projected data. Large variations in va also 
lead to observable effects, thus enabling AMS-02 to be 
sensitive to the details of diffusive reacceleration. 

It is well-known that propagation setups with lower 
Alfven velocities yield lower B/C at energies 1 — 100 
GeV/n [7, 14], and that this can be compensated by an 
increase of L/Dq xx . But, since reacceleration has negli- 
gible influence at high energies, the increase in L/D$ xx 
must be accompanied by a larger value of a so that B/C 
is sufficiently suppressed in the high energy range. On 
the other hand, lower values of va enhance 10 Be/ 9 Be at 
multi-GeV energies which is also compensated by an in- 
crease of L/Dq xx . This behaviour is illustrated in Fig. 4 
where we sketch the 3<r contours from the 10 Be/ 9 Be pro- 
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FIG. 4: The effect of changes in the Alfven velocity, va, in 
the 10 Be/ 9 Be 3<r region projected onto the L vs. Dq xx plane. 
Solid, dashed and dotted lines correspond to va = 36, 15, 
km/s, respectively. Propagation parameters were varied in 

the ranges D 0xx = 0.57 - 19.5 ■ 10 28 cm 2 /s, L = 1.22 - 20.48 
kpc and a = 0.32 — 0.50, and the data set used is that pre- 
sented in Fig. 1. Here, we have neglected convection, and 
marginalized over a. A preference for larger L and lower 
Do xx is evident as va decreases. 

jected data shown in Fig. 1 for different values of va- 
Here, we have used the parameter scan described in the 
previous paragraph, and marginalized over a. 

V. BREAKING THE ASSUMPTIONS 

In order to make the questions addressed in this study 
tractable, we have thus far relied on a number of simplify- 
ing assumptions. In particular, we have assumed homo- 
geneity of the diffusion coefficient D xx over the volume 
of the diffusion zone, considered cylindrical symmetry for 
the system, and adopted a smooth distribution of cosmic 
ray sources with universal injected chemical composition 
and spectra. While such assumptions are reasonable and 
have been useful up to this point in time, they will even- 
tually have to be discarded or revised as they break down 
under the increasing precision of future cosmic ray data. 
In this section, we study a few of the possible ways in 
which data from AMS-02 could potentially require us 
to revise the assumptions commonly made in modeling 
galactic cosmic ray production and propagation. 

We first consider fluctuations in the recent cosmic ray 
injection rate. Cosmic ray sources are indeed believed to 
be of stochastic nature both in space and time — there- 
fore, one does not expect the rate at which cosmic rays 
are introduced into the Milky Way to be homogeneous or 
constant over time. Galactic supernova remnants (SNR), 
in particular, are created at a typical rate of ~ 0.03 per 
year, and stay active for ~ 10 4 — 10 5 yr [33]. Depending 
on whether such an event has taken place recently and 
nearby, the observed cosmic ray spectrum will vary ac- 
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cordingly. The GALPROP code gives the possibility to 
model such stochastic fluctuations [13, 34, 35] by solv- 
ing the transport equation in a three-dimensional spatial 
grid (unlike thus far used in this work) and defining two 
further parameters: the average time tsNR between con- 
secutive SNR events occuring in a kpc 3 volume around 
us, and the time interval ten, during which the SNR 
keeps injecting cosmic rays. We adopt tsNR = 10 4 yr 
— corresponding to ~ 0.03 SNR events per year for a 
standard distribution of sources — , tcR — 10 4 yr and 
run GALPROP with all other parameters as in the true 
model of section III. In the upper plots of Fig. 5 we show 
with thin lines the resulting B/C and 10 Be/ 9 Be for dif- 
ferent positions in the local Galactic Disk. For compar- 
ison, the central thick lines denote the ratios obtained 
at r = 8.5 kpc in the true propagation model with no 
stochastic SNR events (but ran with three spatial dimen- 
sions). Using the thin lines in the upper frames of Fig. 5 
to project AMS-02 data, we show in the lower plots how 
such variations impact the (3cr) propagation parameter 
space inferred. Here, va — 36 km/s, V c $ = dV c /dz = 0, 
the propagation parameters were varied in the ranges 
D 0xx = 4.54 - 8.03 • 10 28 cm 2 /s, L = 3.5 - 6.5 kpc and 
a = 0.39 — 0.43, and we have marginalized over a in the 
lower frames of Fig. 5. Additionally, we provide in Ta- 
ble II the best-fit configurations and reduced chi-squares. 
At some level, stochasticity of cosmic ray sources ulti- 
mately limits our ability to deduce the underlying prop- 
agation model from stable secondary-to-primary and un- 
stable ratio measurements. 

Another example we have considered was the possibil- 
ity that the diffusion coefficient does not follow a simple 
power-law, as assumed in the preceding sections. A sim- 
ple extension of the power-law form is a broken power- 
law, with indices a\ and a>2 below and above a refer- 
ence rigidity, R a . To explore this possibility, we adopt 
a model identical to our previous benchmark model (see 
Sec. Ill), but with diffusion coefficient power-law indices 
of a>i — 0.39 and ot2 = 0.43 below and above the ref- 
erence rigidity, R = {4, 10, 10 2 , 10 3 } GV, while fixing 
D XX (A GV) = 6.04 • 10 28 cm 2 /s. Next, we once again use 
GALPROP to compute cosmic ray fluxes and ratios in 
these models and make projections for the observations 
of AMS-02. Table II shows the results found when we fit 
these projected data using propagation models without 
a broken power-law diffusion coefficient. In particular, 
we use the parameter scan with va = 36 km/s, V Ci o = 
dV c /dz = 0, and ranges D 0xx = 4.54 - 8.03 • 10 28 cm 2 /s, 
L = 3.5 — 6.5 kpc and a = 0.39 — 0.43. Disappoint- 
ingly, good fits are found in each case, thus revealing 
that small variations in a are unlikely to be discern- 
able from AMS-02 data. In order to verify whether 
more extreme variations in a are testable, we consider 
a broken power law for D 0xx with indices a\ — 1/3 and 
a 2 = 1/2, and R = {4, 10, 10 2 , 10 3 } GV, again fixing 
D XX (A GV) = 6.04- 10 28 cm 2 /s. The index 1/3 (1/2) cor- 
responds to a Kolmogorov (Kraichman) power spectrum 
of magnetic inhomogeneities. Proceeding in the same 



fashion as before but using the scan of propagation pa- 
rameters with ranges Dq xx = 0.57 — 19.5 • 10 28 cm 2 /s, 
L = 1.22 - 20.48 kpc and a = 0.32 - 0.50, we obtain the 
best-fit configurations reported in Table II. Acceptable 
fits (and wrong best-fit parameters) result in all cases ex- 
cept for Rq = 10 2 GV. Therefore, if a changes suddenly 
from 1/3 to 1/2 well inside AMS-02 high energy range 
— notice that 10 2 GV corresponds to a proton ( 10 B) ki- 
netic energy per nucleon of 99.1 (49.1) GeV/n — , the 
upcoming data on stable secondary-to-primary and un- 
stable ratios should be sufficient to detect the presence 
of such a drastic break. 

In a similar way, we can test the sensitivity provided 
by the projected AMS-02 data to the source composition 
and distribution. A reasonable possibility is that nearby 
cosmic ray sources produce various species of cosmic rays 
in relative quantities which differ from the average over 
the Milky Way (and differ from the default GALPROP 
assumptions). Unfortunately, Table II shows that we 
have very little sensitivity to a ±20% change in the source 
abundance of carbon nor even to a factor 2 in the source 
abundances of C, N and O. In all cases, very good fits 
to AMS-02 projected B/C and 10 Be/ 9 Be are obtained, 
although the apparent best-fit parameters are not neces- 
sarily the parameters of the true model, opening the pos- 
sibility that we may infer a very well-fit, but incorrect, 
cosmic ray propagation model from AMS-02 data. We il- 
lustrate such effect for the case of a 20% enhanced source 
abundance of C: Fig. 6 shows the 1, 2 and 3cr contours 
using va — 36 km/s, V c .o = dV c /dz — 0, and propagation 
parameters in the ranges D 0xx = 4.54 — 8.03- 10 28 cm 2 /s, 
L = 3.5 — 6.5 kpc and a = 0.39 — 0.43. In each frame we 
have marginalized over the parameter not shown. 

Moreover, we have studied the impact of varying the 
distribution of cosmic ray sources. The usual parame- 
terization optimized to reproduce EGRET gamma ray 
data follows Q inj cx {r/r Q ) ri exp(-^ ! ^ a - j z l pc ) with 
rj = 0.5, £ = 1.0 and a cut-off radius r max = 20 kpc 
[14]. Following Ref. [13], we consider alternative scenar- 
ios, namely a supernova remnant-like distribution (with 
f] = 1.69, £ = 3.33), and a pulsar-like distribution (with 
Q inj cx cosh(r /r c ) exp(- J z 1 | pc )/cosh(r/r c ), r c = 3.5 
kpc). As indicated in Table II, high quality fits are 
found in both cases; and again the best-fit models do 
not coincide with the parameters of the true model. The 
same holds in the case of the usual parameterization with 
rj = 0.5, £ = 1.0 and r max — 20 kpc plus a nearby source 
that we put at (x,y,z) = (8.66,0,0) kpc (Geminga ap- 
proximate position). We are, therefore, forced to con- 
clude that AMS-02 data will likely be insensitive to dif- 
ferent assumptions pertaining to the source composition 
and distribution, and that the values of the propagation 
parameters may even be potentially misinferred as a re- 
sult. 

Now, it is interesting to check if any of the best-fit 
configurations listed in Table II induce cosmic ray fluxes 
that conflict with the corresponding true model. For 
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FIG. 5: The impact of the stochastic nature of sources in the injection of cosmic rays on the B/C and 10 Be/ 9 Be ratios (top 
frames), and on the inferred propagation parameters (bottom frames). In the top frames we show the ratios in different positions 
in the local Galactic Disk — the legend is ordered according to the values of the thin lines at low kinetic energies. In the bottom 
frames we have marginalized over a and the legend indicates the propagation model used to simulate AMS-02 data; in all cases 
va = 36 km/s, V Cj o = dV c /dz — 0, and the propagation parameters were varied in the ranges Dq xx = 4.54 — 8.03 • 10 28 cm 2 /s, 
L — 3.5 — 6.5 kpc and a — 0.39 — 0.43. We have assumed a rate of 0.03 galactic supernova per year and show 3cr contours in 
the lower frames. Stochastic variations limit our ability to deduce the underlying cosmic ray propagation model from stable 
secondary-to-primary and unstable ratio measurements. See the text for more details. 



such, besides B/C and 10 Be/ 9 Be, we also project the 
one-year AMS-02 proton flux, in the range 5 GeV — 1 
TeV and with 3% of assumed systematics. In the last 
column of Table II we show in squared brackets the re- 
duced chi-squares of the combined data set including 
B/C, 10 Be/ 9 Be and the proton flux for the best-fit models 
previously found using B/C and 10 Be/ 9 Be only. In some 
cases — notably when assuming a break from 1/3 to 1 /2 
in the diffusion coefficient index — the proton flux helps 
discriminating wrong non- minimal assumptions. In the 
remaining situations, however, the misinference of prop- 
agation parameters discussed in the previous paragraphs 
still persists. 

Finally, we turn our attention to the source spectral 
index. Up to this point the injection spectrum at the 
sources was assumed to be a double power law in rigid- 



ity with indices "/i = 1.82 and 72 = 2.36, and a break 
i?o = 9 GV as mentioned in section III. We now release 
this assumption by taking several values for 72 — the 
low energy rigidity index 71 and the break Rq are kept 
fixed since we are focussing on high energy data only. 
Proceeding as earlier, we generate the projected AMS- 
02 data for each true model and then find the best-fit 
configuration using only B/C and 10 Be/ 9 Be or including 
also the proton flux — the results are reported in Table 
III, where we have used the parameter scan with ranges 
D 0xx = 0.57 - 19.5 ■ 10 28 cm 2 /s, L = 1.22 - 20.48 kpc 
and a = 0.32 — 0.50. As can be inferred from this Ta- 
ble, the B/C and 10 Be/ 9 Be data set is not sensitive to 
a variation of 72 in the range 2.26—2.46. Such conclu- 
sion was expected because cosmic ray ratios are largely 
independent of the source term. In order to break such 
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B/C + 1U Be/ 9 Be best-fit model (N dof = 21) 
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TABLE II: Best-fit models and reduced chi-squares for the combined B/C+ 10 Be/ 9 Be data set, assuming different variations 
to the underlying true model. For each model, an AMS-02 data set was projected according to the corresponding broken 
assumption. The presented best-fit models were found for va = 36 km/s and V c ,o = dV c /dz = 0. For configurations marked 
with ' + - ) the propagation parameters were varied in the ranges Doxx = 0.57 — 19.5 • 10 28 cm 2 /s, L = 1.22 — 20.48 kpc and 
a — 0.32 — 0.50; for the remaining cases the ranges of the parameter scan were Dq X x = 4.54 — 8.03 • 10 28 cm 2 /s, L = 3.5 — 6.5 
kpc and a = 0.39 — 0.43. Configurations marked with have been run in GALPROP with three spatial dimensions; for 
comparison we also present the best-fit model in the case of generating AMS-02 data with the 3-dimensional true model. In 
the last column the chi-squares in squared brackets were obtained using AMS-02 projected proton flux in addition to the ratios 
B/C and 10 Be/ 9 Be (N dof = 42). 
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TABLE III: Best-fit models and reduced chi-squares for the combined B/C+ 10 Be/ 9 Be and B/C+ 10 Be/ 9 Be-(-p data sets, assum- 
ing different high energy injection indices 72. For each model, the AMS-02 data were projected according to the corresponding 
72. The presented best-fit models were found for va = 36 km/s, V c ,o = dV c /dz = 0, 71 = 1.82, 72 = 2.36 and Ro = 9 GV. The 
propagation parameters were varied in the ranges Dq xx = 0.57 — 19.5 • 10 28 cm 2 /s, L — 1.22 — 20.48 kpc and a = 0.32 — 0.50. 



degenaracy a cosmic ray flux measurement — e.g. of pro- 
tons — must be used. The proton flux at high energies 
is supposed to go as _E~( Q +T 2 ) anc l consequently it helps 
discriminating wrong values for the injection index 72. 
In fact, as shown in Table III, the analysis with a com- 
bined data set including B/C, 10 Be/ 9 Be and the proton 
flux produces relatively large reduced chi-squares of the 
best-fit configurations. This indicates that 72 is devi- 
ating from the assumed value of 2.36. Notice as well 
that, when one uses the projected proton data, the best- 



fit values for a are smaller for smaller injection indices 
72, and vice- versa — this is because the proton flux is 
sensitive to the combination a + 72. To proceed further 
and quantify how sensitive AMS-02 will be to the cos- 
mic ray injection spectrum, we scan the parameter space 
(D axx ,j 2 ,a) and fix v A = 36 km/s, V Cj0 = dV c /dz = 0, 
and L = 5 kpc (notice that D 0xx and L are approxi- 
mately degenerate). We ran GALPROP 343 times, in a 
7x7x7 grid of the parameters (Dq xx , 72, a) over the ranges 
D 0xx = 4.54 - 8.03 • 10 28 cm 2 /s, 72 = 2.33 - 2.39 and 
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FIG. 6: Regions consistent (within 1, 2 and 3er) with pro- 
jected B/C and 10 Be/ 9 Be AMS-02 data in the L vs. D 0xx , 
a vs. Dqxx, and a vs. L planes. Here, the AMS-02 data 
were projected assuming a true model with a source abun- 
dance of 12 C 20% higher than in the default model of Ref. [6]. 
Propagation parameters were varied in the ranges Dq X x = 
4.54 - 8.03- 10 28 cm 2 /s, L = 3.5-6.5 kpc and a = 0.39-0.43. 
We have assumed va = 36 km/s, V c ,o = dV c /dz — 0, and have 
marginalized in each frame over the parameter not shown. 
Notice that, although a good fit was found to the data in 
this case, the apparent best-fit parameters are significantly 
different from the values of the true model. 



FIG. 7: Inferred parameter regions (within 3<r) in the plane a 
vs. 72. The dashed (dotted) [solid] lines refer to the case where 
the AMS-02 projected data set consisted of the proton flux 
(B/C and 10 Be/ 9 Be) [proton flux, B/C and 10 Be/ 9 Be]. Here, 
the propagation parameters were varied in the ranges Dq X x = 
4.54 - 8.03 • 10 28 cm 2 /s, 72 = 2.33 - 2.39 and a = 0.38 - 0.44. 
We have assumed va = 36 km/s, V c ,o ~ dV c /dz — 0, and 
have marginalized over Do xx - The complementarity between 
the ratios B/C and 10 Be/ 9 Be and the proton flux is evident. 

a = 0.38 — 0.44. Linearly-(Logarithmically-)spaced grid- 
points were implemented for 72, a (Dq xx ). An infill of 
3 points between consecutive gridpoints was applied re- 
sulting in 25x25x25=15,625 different models. Figure 7 
shows the 3tr regions in the a vs. 72 plane (marginalized 
over Dq xx ) from the projected AMS-02 measurements of 
the (i) proton flux (dashed lines), (ii) B/C and 10 Be/ 9 Be 
(dotted lines), and (Hi) B/C, 10 Be/ 9 Be and proton flux 
(solid lines). Such as anticipated earlier, the ratios B/C 
and 10 Be/ 9 Be are efficient probes of the diffusion index 
but insensitive to 72, while the proton flux constrains 
essentially the quantity a + 72. Hence, within minimal 
propagation models, AMS-02 has the potential to pin- 
point both the diffusion index a and the high energy 
source spectral index 72 with good accuracy. Note that 
this result stems exactly from the combination of the B/C 
and 10 Be/ 9 Be ratios and the proton flux. While a reason- 
ably precise high energy proton flux has been measured 
by past and present instruments (e.g. PAMELA [36]) — 
allowing corrispondingly precise estimates of a + 72 — , 
only AMS-02 (or future detectors) will be able to break 
down the degeneracy between a and 72 with high energy 
quality cosmic ray ratios. 

VI. CONCLUSIONS 

In this paper, we have considered the ability of the up- 
coming AMS-02 experiment to measure stable secondary- 
to-primary and unstable ratios (such as boron-to-carbon 
and beryllium-10-to-beryllium-9), and studied to what 
extent this information could be used to constrain the 
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model describing cosmic ray propagation in the Milky 
Way. Within the context of relatively simple propaga- 
tion models, we find that the parameters can be very 
tightly constrained by the projected AMS-02 data; con- 
siderably more so than is possible with currently existing 
data [6, 7]. 

The ability of AMS-02 to constrain the cosmic ray 
propagation model can be considerably reduced, how- 
ever, if more complex models with larger numbers of free 
parameters are considered. Only a rough determination 
can be made of the Alfven velocity (which dictates diffu- 
sive reacceleration) in most cases, for example. On the 
other hand, if even a relatively small degree of convection 
(at the level of ~10 km/s/kpc) is present, this is expected 
to be discernable from the AMS-02 data. Other aspects 
of cosmic ray propagation (including, for example, the 
detailed energy dependence of the diffusion coefficient, or 
variations in the source distribution or injected chemical 
composition) are unlikely to be significantly constrained 
by upcoming data on stable secondary-to-primary and 
unstable ratios. In some cases, we have found that the 
parameters of the underlying diffusion model could be 
misinferred due to inaccurate assumptions implicit in 
the model. Local variations in recent supernova activ- 
ity could also lead to somewhat skewed determinations 
of propagation parameters. The source spectral index, in 
constrast, is likely to be well constrained within minimal 



models if proton flux measurements are used in addition 
to the ratios B/C and 10 Be/ 9 Be. With the latter case we 
have therefore exemplified how in some situations cosmic 
ray observables other than ratios are powerful tools in 
discriminating non-minimal assumptions. 

In summary, the introduction of data from AMS-02 
will make it possible to significantly expand our under- 
standing of how cosmic rays propagate through the in- 
terstellar medium of the Milky Way. As we have demon- 
strated, the characteristics of simple propagation models 
will be tightly constrained by this data set. Moving be- 
yond such simple scenarios, some of the underlying model 
assumptions will also become testable with such data, al- 
lowing us to better refine our cosmic ray predictions, and 
more generally expand our understanding of the Galactic 
cosmic ray spectrum. 
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